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In the limit when the two-body scattering length a is neg- 
ative and much larger than the effective two-body interac- 
tion radius the contribution to the ground state energy due 
to the three-body correlations is given by the Efimov effect. 
For particular values of the diluteness parameter p|a|^ the 
three-body contribution can become the dominant term of 
the energy density functional. Under these conditions both 
Bose-Einstein and Fermi-Dirac systems could become self- 
bound and either boson droplets or fermion "designer nuclei" 
of various sizes and densities could be manufactured. 

PACS numbers: 03.75.Fi, 05.30.Jp 

The extraordinary progress in realizing Bose-Einstein 
Condensates (EEC) and of the Fermi-Dirac counter parts 
in traps |^,^ raises the question of whether they can exist 
also in isolation, without a trap potential. In the dilute 
limit, the physical properties are entirely determined by 
the two-particle scattering length a. To leading order in 
a the energy per particle of a uniform condensate is given 
by 
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where p is the number density and m is the particle mass. 
For positive a the energy is positive, suggesting that there 
is no binding. The energy density, including higher order 
corrections, is given by 
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where C is a constant P-H|. These higher order correc- 
tions depend on dimensionless parameter pa^ up to the 
term C, which depends on the range of the two-particle 
interaction. A similar expression to Eq. (y), but for the 
energy of a dilute Fermi system, in the case of positive 
scattering lengths, was derived essentially in parallel with 
the Bose case [y,^]. The positive energy state is stable 
with respect to long wave density fluctuations, but is only 
metastable if the underlying two-particle interaction sup- 
ports bound states. In fact most of the BEC decay into 
the two-particle bound states, but they are still quite 
long lived. The decay rate of the condensate depends 
on a as r oc a'*p^, see Refs. M and earlier references 



cited therein for details and some rather unexpected fea- 
tures. In contrast, when the effective scattering length 
is attractive (a < 0), the energy expression Eq. p) is 
unbounded from below and the uniform density state is 
unstable with respect to density fluctuations. The time 
scale here is very short (F ex ■\/|a|p). It would thus seem 
that a droplet of BE condensed matter would be either 
unbound or would decay very quickly. This is not nec- 
essarily the case. The picture is radically changed when 
one considers the effects of the three-body correlations. 
Recently it was confirmed experimentally the fact that 
one can now control the two-body scattering length by 
placing the atoms in an external field 0], see also Refs. 
[Q. For those cases where successful Bose-Einstein con- 
densation was achieved, the bare or unadulterated by any 
external fields two-body scattering length a was compa- 
rable in magnitude with the two-body potential radius. 
More exactly, the value of the scattering length is compa- 
rable with the van der Waals length [gj^ . A new regime 
is now attainable, when the Bose system is still dilute, 
but with respect to a new characteristic length. If in 
the two-body system the scattering length is very large, 
i.e. if \a\ 3> rg, there are now two independent dimen- 
sionless parameters, pa^ and pr^. In this limit in the 
three-body system one has a very unusual phenomenon, 
the so called Efimov effect fiq], which is manifested in 
the appearance of a very large number of three-body 
bound states N « sq/tt ln(|a|/rQ), where sq « 1.0064, 
all with the same quantum numbers 0^. The spatial ex- 
tensions of these Efimov states range from distances of 
the order of 0(ro), for the tightest bound one, to dis- 
tances 0{a), for the least bound state. These three- 
body bound states appear irrespective of whether the 
two-body scattering length is positive or negative. The 
sizes of these states change from one state to the next by 
a factor exp(7r/so) ~ 22.68, while their energies change 
by a factor exp(— 27r/so). The spectrum has an exponen- 
tial character and the bound state wave functions obey a 
simple scaling law. The properties of these states easily 
follow from the fact that in the region rg <C R <C |a|, 
where R^ — 2(r^2 + r^^ + ''3i)/3 and r^ is the distance 
between particles k and I, (k,l — 1,2,3), there is an 
effective three-particle interaction with an universal at- 
tractive character —SpTi /2mR^. The universality of this 
effect resides in the fact that all its properties are fully 
determined by the scattering length a, the radius of the 
two-body interaction ro, the logarithmic derivative A of 
the three-body wave function at i? « rg and the dimen- 
sionless constant Sg. If the particles have spin, the sit- 



uation can become a little bit more complicated, as the 
Efimov states can appear for various values of the total 
spin of the three-body system. Typically the spectrum 
ceases to be strictly exponential and the wave functions 
have a somewhat more complex structure. However, on 
average the qualitative features of these new type of Efi- 
mov three-body states remain largely unchanged |0] . 

Efimov |13| has also shown that at zero energy the 
amplitude for the three-particle collisions is determined 
by the same universal interaction —s^fi /2mR^. In the 
case of negative scattering lengths (a < 0) this amplitude 
is given by the following universal formula 
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where the numerical values of the universal constants di , 
d2 < and oq have been determined numerically recently 
by Braaten et al. |0] . Here oq is the value of the two- 
body scattering length for which the first three-body 
bound state with zero energy is formed. In the present 
paramctrization oq replaces the logarithmic derivative 
and the matching radius used by Efimov |12,nq]. Un- 
like di and d2, the parameter ao is system dependent 
and is also a genuine three-body characteristic. When 
approaching the three-body threshold (by changing the 
strength of the two-body interaction), just before the 
three-body bound state is formed, g^ — > — oo, in com- 
plete analogy with the behavior of a two-body scatter- 
ing amplitude. After the three-body state has appeared 
(73 is positive and 33 — > oo when the threshold for the 
appearance of the three-body bound state is approached 
from the other side. 

If the three-body zero-energy scattering amplitude 33 
is known, the contribution of the three-particle collisions 
to the ground state energy density of a dilute Bose gas 
can be evaluated in a similar manner as the leading order 
contribution of the two-body collisions 
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The experience gained from studying corrections to the 
main two-body contribution to the ground state energy 
of a dilute quantum gas, see Refs. |^,H,Q and Eq. (g), 
shows that all other contributions are small in the dilute 
limit and they are not expected to lead to any qualitative 
changes in the properties of these systems. (In the case of 
fermion systems pairing correlations however could lead 
to a significant rearrangement of the ground state prop- 
erties, but energetically the correction is typically small 
and the density profiles are not modified in any drastic 
manner, see e.g. nuclei uM.) The Szip) contribution to 
the energy density can dominate the contribution aris- 
ing from the two-body collisions if the argument of the 
tangent is close to {2n+ l)7r/2 and the scattering length 
a « ao exp(n7r/so), n = 0, 1, ... is such that a three- 
body state is on the verge of appearing or it has just been 



formed. The situation is somewhat unique in this limit. 
With respect to the two-body collisions the system is ex- 
tremely dilute, but somewhat less dilute with respect to 
three-body collisions. 

At the points where a three-body bound state appears, 
where (73 becomes infinite and in the immediate neigh- 
borhood of them, the present calculational scheme fails, 
since the contribution of the three-body collisions has 
been evaluated only in the leading order of the gas ap- 
proximation with respect to the three-body collisions. 
In the region near such poles, the contribution of the 
three-particle collisions could dominate over the genuine 
two-particle contributions, for appropriate values of p 
and a. In a trap the density profile of such a trapped 
Bose-Einstein condensed gas is given now by the new 
"Thomas-Fermi" formula 
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where Vext {f) is the trapping potential and p is the chem- 
ical potential. Besides obvious changes in the density 
profile of a trapped Bose-Einstein condensed gas, the 
spectrum of the elementary and collective excitations is, 
naturally, modified as well, as the compressibility of such 
a system is significantly affected by the three-body con- 
tribution to the ground state energy £3. I would like 
to note here that the analyses of Refs. |15[ consider a 
formally similar situation, but with a fictitious repulsive 
three-body force, whose nature and strength are never 
specified. 

Since the two-body contribution to the ground state 
energy of a dilute Bose gas is negative, the three-body 
collisions in the regime where 33 > could lead to the 
stabilization of the system. What is particularly interest- 
ing for such a system is that a boson droplet - a boselet 
- could become self-bound and the trapping potential is 
not required anymore to keep the particles together. In 
the absence of the trapping potential and for a very large 
number of particles a boselet will have an almost constant 
density, corresponding to the infinite matter equilibrium 
density 
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The ground state energy of an ensemble of N bosons in 
the meanfield approximation is given by |16|] 



EiN) = / dVe(r) = / dV 



2m 



:92p{r) 



6 



93p{rf 



(7) 



where e{r) is the energy density and p{r) — |f/'(^)P is 
the number density. The density profile of semi-infinite 
matter can be shown to be given by 



p{z) = \^{zt 



Po 



1 + cxp(2ko-2) ' 



where 



Kg 



'2m|/io| 



Mo 



M 
833' 



(8) 



(9) 



z is the spatial coordinate normal to the surface and /io 
is the chemical potential corresponding to infinite matter 
at equilibrium density po- One can now determine the 
surface tension a of boselets from the obvious relation 
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The spectrum of both volume and surface sound waves 
can then be easily specified. The density profile of an 
infinite slab of finite width has a simple expression as 
well 
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where p = —h k'^ /2'm and the chemical potential satisfies 
the restrictions mo < M < 0- It is straightforward to 
show that for any slab with a width larger than its skin 
thickness, the quantities (m — Mo)/mo £ind (kq — k,)/kq 
are both exponentially small. These facts (along with 
numerical evidence not presented here) suggest that the 
basic properties of a A^-particle boselet are given by the 
following relations 



p{r) « Mo 1 + 



1 



2kR 



1 



cosh(2Kr) ' 



cosh(2Ki?) 

-1/3 _^ 0(^-2/3)^ 

E{N) = fioN + 47rrgcrA^2/3 _^ oi^^Us^^ 



R = roN 



1/3 



riN 



al{l-l){l + 2) 
mpoR? 



+ 0(Ar-4/3)^ 



(12) 

(13) 
(14) 

(15) 



where Tq = (3/47^/9o)^/^ M = dE{N)/dN = ~r?K^l2m 
and LDi is the frequency of the surface vibration mode 
with angular momentum I. The absence of the constant 
term in the expression for the radius was established nu- 
merically. The central density is larger then po due to 
surface tension. 

The possibility that the entire system can also undergo 
a transition to a gas phase of trimers cannot be ruled out 
at this time. Since there are no Efimov states for four or 
more particles this trimer phase is perhaps unique. One 
cannot fail to see here an analogy with the Cooper pair- 
BEC crossover in the fermion case |1^. At each new 
three-body threshold, when 93 becomes infinite, a new 
trimer phase appears, made of spatially larger trimers. 
The density drops naturally by a factor of 3 if a trimer 



phase is formed. The trimer-trimer scattering length is 
expected to be of the order of the trimer size, i.e. of 
order a, but the sign of this trimer-trimer scattering am- 
plitude is so far unknown. Upon collapsing into trimers 
the interaction energy decreases significantly, as now this 
energy is controlled by effective two-body processes only. 
If a trimer phase is formed, the size of the cloud in the 
trap should change abruptly. 

It is obvious that in the case of a Fermi-Dirac system 
the role of three-body collisions could play an analogous 
role, if the Efimov effect could take place. Since three 
identical fermions could never be all in a relative s-state, 
the particles should have either a spin larger than 1/2 
or some other additional degree of freedom, like isospin 
in the case of nucleons. The characteristics of the Efi- 
mov effect for the case of particles with arbitrary spin 
and/or other discrete degrees of freedom (as well as arbi- 
trary masses) have been established in Ref. [Q . Let me 
consider for illustrative purposes the case of two identi- 
cal spin-1/2 fermion species tt and v. In the meanfield 
approximation the ground state energy can be evaluated 
using the following energy density functional 

K(r)+M.(r)]^-|[p.(r)2+p,(r)2] 
p^ir)p,{r)[p^ir) + p,{r)]Y (16) 
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where p-n^vi'i') are the corresponding number density 
distributions, Nt^^^ — ^ d^rpT^^^{r) and ?i^r^,i/(r) = 

^^ Sn='i 2|VV'n '"'(»") P are the momentum distribu- 
tions, defined through the corresponding single-particle 

wave functions tpn'^ii")- 

Again, the most interesting regime to consider is that 
of a negative two-body scattering length 32 < 0. If the 
three-body term is repulsive, that is if 173 > 0, a dilute 
Fermi-Dirac droplet will behave very much like a nu- 
clear system. In this case the two-body effects will have 
an attractive role and the energy density functional of 
such a system will have the same qualitative structure as 
the popular Skyrme energy density functional in nuclear 
physics ||l^. The Fermi-Dirac droplets - the fermilets 
- will have entirely unexpected properties: they will be 
self-bound and show saturation properties as well, see 
Fig. 1. The existence of an equilibrium state for infi- 
nite homogeneous fermionic matter is a sufficient condi- 
tion for the existence of finite systems as well, see the 
case of boselets for example. Since one would be able 
to change the relative magnitude of the two-body and 
three-body interactions the central density of a fermilet 
can be controlled. In the absence of Coulomb interaction 
the number of particles in a single fermilet is arbitrary. 
A particularly interesting aspect could be the interplay 
between the formation of Cooper pairs and fermionic 



trimers, since when 173 > a trimer bound state exists. 
The possibiUty to have a Fermionic system, with attrac- 
tive two-body effective interactions and repulsive three- 
body effective interactions, opens thus the way towards 
the creation of "designer nuclei" , an almost unthinkable 
flexibility, which could not be matched even by atomic 
clusters. Self-bound droplets of mixtures of fcrmions and 
bosons - ferbolets - are expected as well. 

As a final remark, since recombination now requires 
four particle to collide, it is not unreasonable to expect 
rather long lifetimes for these new objects. Three-body 
recombination into deep two-body bound states (if such 
states exist) could however define the lifetime of these 
objects B. 
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FIG. 1. A typical contour plot of the energy density for an 
homogeneous Fermi system consisting of two fermion species 
TT and u, see Eq.jlq), for the case h — m — gs = 1 and 
g2 — —5. Only the negative part of the energy density surface 
is plotted. The local minimum and the two saddle points are 
labeled by M and S respectively. For other sets of parameters 
there is either only one saddle point or none at all. 
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